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Abstract: Using lattice simulations, we measure the sphaleron rate in the Standard Model 
as a function of temperature through the electroweak cross-over, for the Higgs masses 
rriH = 115 and mn = 160 GeV. We pay special attention to the shutting off of the baryon 
rate as the temperature is lowered. This quantity enters computations of Baryogenesis 
via Leptogenesis, where non-zero lepton number is converted into non-zero baryon number 
by equilibrium sphaleron transitions. Combining existing numerical methods applicable 
in the symmetric and broken electroweak phases, we find the temperature dependence of 
the sphaleron rate at very high temperature, through the electroweak cross-over transition, 
and deep into the broken phase. 
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1 Introduction 

As a result of a quantum anomaly, baryon number B and the lepton numbers Lj are not 
strictly conserved charges in the Standard Model. As the Chern-Simons number Acs of 
the gauge field evolves in time, because only left-handed fermions are coupled to the SU(2) 
gauge field, the baryon and lepton numbers B and Lj change according to the relation (see 
for instance [1, 2]) 

— [B{t) - B{0)] = Li{t) - L,(0) = Acs(t) - Acs(O) = f dt [ x Tt F^" F^, , 
no 327r^ Jo J 

(1.1) 

where F^'^ is the field strength tensor of the SU(2) gauge field and no = 3 is the number 
of fermion generations in the Standard Model. 

In addition, the electroweak sector of the Standard Model has an infinite set of "pure 
gauge" degenerate vacua, corresponding to integer values Acs- The question of baryon and 
lepton number violation in the Standard Model therefore becomes a question of whether 
dynamics allow transitions from one vacuum to another. 
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At zero temperature transitions occur via instantons [3], and the rate is minute, neg- 
ligible even on cosmological scales. However at finite temperature, thermal fluctuations 
can lead to spontaneous transitions, controlled in equilibrium by a diffusion (or sphaleron) 
rate Tam, 

V, t-^oo V t 

At zero temperature, the Standard Model Higgs field acquires a vacuum expectation value 
= v/^/2, with V = 246 GeV. In contrast, at high temperature {\(f)\) ~ 0, and there is 
a "symmetry breaking" transition between the two phases ("symmetric" and "broken"). 
Note that, because of gauge symmetry and Elitzur's theorem, (i;^) is always zero, whereas 
(101) is always non-zero (the value depends on the renormalization scheme, though). Thus, 
we do not have an exact local order parameter which would vanish in the symmetric phase 
and be non- vanishing in the broken phase. Indeed, at low Higgs masses (m/f < 72 GeV) , 
the Standard Model phase transition is an actual first-order phase transition, but at exper- 
imentally allowed values for the mass the electroweak transition in the Minimal Standard 
Model is a continuous cross-over [4-7]. Nevertheless, we shall use the labels "symmetric" 
and "broken" phases to describe the states above and below the cross-over temperature. 

In the symmetric phase the sphaleron rate is unsuppressed and is proportional to 
T'^. In the broken phase the energy barrier between the topological sectors grows as 
increases, and parametrically the rate is 

— — - i^hxoVOL^e. T , [1.6) 

where £'sph is the sphaleron energy (energy barrier) and Kbrok is a numerical coefficient [8] . 
In the physical cross-over case it is of interest to know how the rate "shuts off" as tem- 
perature decreases, including the precise determination of the rate in the exponentially 
suppressed temperature range. This is important e.g. for Leptogenesis scenarios, as dis- 
cussed in Sect. 5. 

1.1 Calculations of the sphaleron rate 

The sphaleron rate has been computed extensively using lattice simulations. For pure 
SU(2) gauge theory, the quantum diffusion rate is approximately recovered in the classical 
dynamics of a dimensionally reduced effective theory [9-12], taking into account the proper 
conversion factors [13, 14]. The calculation was improved by including the effect of hard 
thermal loops [15, 16], and the magnitude was finally settled using the Langevin dynamics 
of another effective theory [17, 18], again using a conversion prescription [18-20]. 

The rate in the electroweak theory (i.e. including the Higgs scalar) was computed in 
the symmetric phase at mn — 34 GeV in [11], and for mn ~ mw in [21] in approximate 
agreement with the pure-gauge result. The most recent simulations [22] combine all previ- 
ous methods, extends the effective theory of [17] to include the Higgs field, and finds that 
at niH ^ 44 GeV [22], 

!^.(8..4±„.10)(^)(ln(p|)+c).L (1.4, 
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(1.5) 



with ^ 

ml = ^9^T\ C = 3.041, = f-. 

D 47r 

In the broken phase, the rate is less well known. The straightforward numerical methods 

of [11, 21] are not able to cope with the large suppression of the rate, a problem resolved by 

Moore using multicanonical methods [23, 24]. In these papers the low Higgs mass region, 

where the phase transition is strongly first order, was considered. The latter paper finds 

that at the transition temperature Tc the broken phase rate is 



"pbrok/'T^ 
diff \-'^ 



'-pA 
c 



exp 



-R 



GeV/ 



(1.6) 



with 



i2(49.2) = 24.7 ± 0.4, i?(44.8) = 28.3 ± 0.4, i?(41.2) = 31.2 ± 0.6. (1.7) 



In this case the strong suppression in the broken phase results in an effectively instantaneous 
shut-off of sphaleron processes during the strong first-order phase transition. In the cross- 
over regime, one would expect a gradual shut-off as the Higgs field expectation value 
increases. 

In [25] the sphaleron rate was calculated in the electroweak cross-over region, around 
rriH = 120 GeV, from high-temperature into the low-temperature phase, using classical 
dynamics. Shortly afterwards [22], a similar calculation was performed at mn = 130 GeV, 
but using the effective dynamics of [17]. In both cases a rather rapid shut-off of sphaleron 
diffusion through the cross-over was reported. However, because no multicanonical methods 
were used, the exponentially suppressed tail was not resolved in detail. The aim of this 
paper is to improve on these calculations, and pin down the sphaleron rate through the 
cross-over, also in the domain of exponential suppression. 

Some time ago [26], known calculations of the sphaleron rate were collected and ex- 
trapolated to the cross-over region. The conclusion was that within a range 



and for a range of T so that 



we have 



-In 



2^4 



100 GeV <mH< 200 GeV, 
-In [rdifr(r)/r4] ~ 30- 50, 

(iTiH - 150 GeV V /T- ISOGeV^-'' 



10 GeV 



Y 



10 GeV / 



(1.8) 



(1.9) 



with Coo = 39.6, cio = 3.52, cqi = -7.09, C20 = -0.376, cn = 0.421, C02 = 0.17. We will 
use this result for guidance and establish its range of validity. 

In this work we study the sphaleron rate at Higgs mass values mn = 115 GeV and 
rriH = 160 GeV^, with special attention paid at the rate deep in the broken phase. 



^ These Higgs mass values were within the experimentally allowed range when the simulations were 
started [27]. 
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1.2 The sphaleron rate and Baryogenesis 



In Electroweak Baryogenesis [1, 2], bubble nucleation in a first-order electroweak phase 
transition divides space into regions of different electroweak phases. Inside the bubbles 
is the broken phase where the sphaleron rate is very small. Outside the bubbles is the 
symmetric phase, where the rate is large. CP-violating interactions between the fermions 
in the plasma and the advancing bubble walls generate a net chiral fermion current, which 
is then transformed, by sphaleron processes, into a baryon and lepton asymmetry [28]. 
A first-order phase transition is not realized in the Minimal Standard Model, but the 
Baryogenesis scenario may be relevant for extensions such as the 2-IIiggs Doublet Model 
and Supersymmetric Standard Models, in which case a reliable calculation of the sphaleron 
rate in both phases is important. 

In Baryogenesis via Leptogenesis [29, 30], a lepton asymmetry is assumed to originate 
from some separate process, represented here by a time-dependent source fi{t) which may or 
may not be active at the electroweak scale. Sphaleron transitions equilibrate the system, so 
that the lepton asymmetry is transformed into net Li and B. Following [26], the equations 
controlling this equilibration read 



where T = T{t), ^{t) = 7[rdifr(r), 'y(r)], r](t) = r][v{T)] and no = 3. v{T) is the expecta- 
tion value of the Higgs field, taken to be v{T) ~ (0^)- 

In this paper, we will calculate rdiff(T) and v{T), in the Minimal Standard Model at 
experimentally allowed Higgs masses, i.e. in the regime where the electroweak transition is 
an equilibrium cross-over. 

In Section 2 we briefiy explain how to treat the dynamics of the Standard Model 
at finite temperature as an effective 3-dimensional SU(2)-Higgs theory. Section 3 we de- 
scribe the numerical lattice Monte-Carlo methods employed here, in the high- and low- 
temperature regimes, respectively. In Section 4 we present our results for rdifr(T) and 
v{T). In Section 5 we solve (5.2) for some simple cases to assess the impact of including 
the correct sphaleron rate rather than assuming an instantaneous shut-off. We conclude in 
Section 6. 

2 SU(2)-Higgs model on the lattice 

The sphalerons are non-perturbative field configurations and thus we need to use non- 
perturbative lattice simulations in order to calculate the rate reliably. Furthermore, because 
of the infrared problem in the thermodynamics of the non-Abelian Yang-Mills theory, the 
modes with momenta k < g^T are non-perturbative. This means that the perturbative 
expansion becomes impossible beyond some loop order [31]. 




(1.10) 



(1.11) 
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2.1 Dimensional reduction 

The full four-dimensional Standard Model with chirally coupled fermions is too unwieldy 
to simulate on the lattice. However, for static (time-independent) thermodynamics we 
can use the fact that the weak coupling constant is small, and apply perturbation theory 
only to modes which can be reliably treated with perturbative methods: that is, to modes 
with momentum k > g^T, where is the weak gauge coupling. This procedure is called 
dimensional reduction, because it results in a three-dimensional effective theory for the 
soft (g'^T) modes. The effective theory is purely bosonic, and it fully includes the essential 
non-perturbative physics. The detailed description how this is performed can be found in 
ref. [32-34]; for earlier and related work, see [35-40]. 

The perturbative derivation of the effective theory is based on the hierarchy between 
the hard {k>T), electric {k ~ gT) and magnetic {k ~ g^T) scales on an Euclidean finite- 
temperature path integral. In the first stage we integrate over the hard scales, obtaining an 
effective theory of scales k < gT. Because the Matsubara frequencies for the bosonic and 
fermionic field modes are kQ°^°^ = 27rnT and k^^^^°^ = 7r(2n -|- 1)T, n ^ Z all fermionic 
modes and non-static (/cq 0) bosonic modes are of order T. Thus, the effective theory 
is purely bosonic and three-dimensional. Concretely, the actual "integration" is done by 
writing down a general renormalizable effective theory and matching the perturbatively 
computed two-, three- and four-point functions in the effective theory and in the original 
four-dimensional theory, thus fixing the parameters of the effective theory. 

The effective theory can be further simplified by integrating over scales gT, which gives 
us the three-dimensional SU(2) gauge theory coupled to a scalar (Higgs) field:^ 

S= j d?x Qi^^i^:;- + {Dick)^ (D.cP) + mlcpU + WUf^ ■ (2.1) 

The coupling constants of the theory are 53 ~ g'^T, and A3. These depend on tempera- 
ture and the parameters of the Standard Model; the full expressions are given in ref. [32]. It 
is customary to use dimensionless quantities x and y and express the set of the parameters 
of the effective theory as 

9l x = ^, y=4- (2.2) 

Here the dimensionful parameter g'^ gives the scale and the physics is completely determined 
by the values of x and y. 

The values of x and y for the Higgs masses used in this work, mn = 115 GeV and 
m-H = 160 GeV, are shown in Fig. 1 over the relevant temperature range. The other 
significant Standard Model parameters which influence the values of x and y are the Fermi 
coupling Gf = 1.16639 x 10-^ mz = 91.1876 GeV, mw = 80.425 GeV, nit = 174.3 GeV 
and as{mw) = 0.1187. 



^The hypercharge U(l) and gluon SU(3) gauge fields in principle survive the dimensional reduction. 
However, the gluons do not couple to the weak gauge and Higgs fields and only modify the parameters in 
Eq. (2.1) through loop corrections. The U(l) field is omitted here because its effect on the transition is 
numerically small [4]. 
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Figure 1. The values of a; (left) and y (right) for ruH = 115 and 160 GcV in the temperature range 
of interest. 



2.2 Lattice-continuum relations 

At this point we have built a continuum 3D theory from the 4D fundamental one. Now we 
have to put it on the lattice. We use here the common lattice discretization of the action [5] 



+ ^ Tr[$t(x)$(x)j +f3nY 



-Tr 



$t(x)$(2;) 



(2.3) 



which is constructed only with gauge-invariant terms. Here Ui{x) is the SU(2) gauge link 
variable, Pij is the standard ij-plane plaquette constructed from the link variables, and the 
lattice Higgs field is naively scaled from the continuum field with (l/8)/3G/3_f/$^^' = 
The parameters /3g, Ph and /Sr are related to the parameters g^a, x and y by 



9l 



gia 



/32 



H 



8 [(3h 



+ — (1 + 4x) + 



Pg J 



327r 



167r2 



51 
16 



+ 9x- 12x' 



In^ + C ) +5.0 + 5.2X 



(2.4) 
(2.5) 

(2.6) 



Here S = 3.1759115 and C = 0.08849 and the other numerical constants appearing in (2.6) 
are specific for the SU(2) + Higgs model and are computed in [41]. 

We note here that because the action is superrenormalizable, only the mass term (y) 
gets renormalized. The continuum limit is taken using equations (2.4)-(2.6) by letting 



/3g — ^ oo while keeping x, y and constant. The counterterms in Eq. (2.6) remove the 
Unear and logarithmic divergences in the lattice spacing a. 

We actually implement the following simple improvement to equations (2.4)-(2.6), 
which cancels part of the 0{a) corrections [42]: in the lattice action (2.3) and equa- 
tions (2.5)-(2.6) we substitute 

/3g ^ /3G,improved = ^/{glo) + 0.6674 , (2.7) 

which makes the gauge part of the action 0{a) accurate. In fact, (2.5) can also be 0{a) 
improved [42], but since the full 0{a) correction to (2.6) is not known we do not implement 
it here. In what follows we consistently use the identity /3g = 4/ ((73 a). 

The gauge invariant Higgs field expectation value ($^$) has linear and logarithmic 
UV divergences. Subtracting the divergences from the lattice expectation value, we obtain 
the continuum quantity ((/>^(/>) as [5] using 

. ((.t,, - ^) - ^ (,o. ^ . , , i.. _ ,) . O (i-) ,.B, 



where C, + \T? - 6 ^ 0.6678. 
2.3 Real-time evolution 

The effective theory in Eq. (2.1) is well understood and has been very successfully used in 
studies of static thermodynamical quantities of hot electroweak physics. As such, it does 
not describe dynamical phenomena, which include sphaleron transitions. It is possible to 
take the theory in Eq. (2.1) and use the classical equations of motion to describe the time 
evolution of the fields, as was done in the early studies of the sphaleron rate [11, 21, 23]. 
However, it has been shown that the classical theory contains divergent UV contributions 
to the gauge field dynamics, and the results are cut-off dependent [43]. Hence, technically 
the infrared gauge field dynamics of the classical theory do not exist. The physical origin 
of the problem is that the Landau damping of the transverse gauge fields in the classical 
theory is UV divergent, and the theory does not have a physical continuum limit. 

These problems can be ameliorated by studying classical theory with hard thermal loop 
(HTL) effects included [44]. This leads to complicated and expensive numerical implemen- 
tations [15, 16]. However, as first demonstrated by Bodeker [17], the physical damping 
makes the dynamics of the infrared gauge field modes (modes with k > g^T) to be fully 
overdamped. Then, at leading order in l/\Q.{l/g) the evolution of these modes is described 
with simple Langevin dynamics (in = gauge, and identifying H/T = S in (2.1)) [20]: 



r) J-f 

dtA. = -a;,'— + Ct, (2.9) 



where a^i is the non-Abelian "color" conductivity. 



a- = % (,„!=£ + 3^041), (2.10) 

mj) 47r \ 7 ' 
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and again = (11/6)(7^T^ in the standard model. ^ is a random Gaussian noise with 

(Cf (x,t)e}(x',t')) = 2a,iT6,,5^'6i^ - x')<5(t - t'). (2.11) 

The Higgs field has parametricahy much less damping. Hence, it can also be evolved with 
a Langevin equation, but with a much faster rate of evolution. To this accuracy we can 
take it to be infinitely fast in comparison with the gauge field evolution [22]. Iterating, we 
can solve for 7 = 0.66361688 g^T. 

In principle the Langevin evolution is straightforward to implement on the lattice. 
However, it is unnecessarily slow: we can substitute it with any dissipative update, as long 
as the relation between the evolution rates is known. Thus, it is much more efficient to use 
random-order heat-bath update algorithm for the SU(2) gauge fields [18, 45]. Now n full 
heat-bath update sweeps through the lattice correspond to the real-time step 

At = ^ac\a^n. (2.12) 

We note that this relation is valid for "unimproved" single-plaquette Wilson gauge and 
random-order heat bath; for other choices the relation would differ. The Higgs field is 
updated with a mixture of heat bath and overrelaxation much more frequently than the 
gauge field [45]. 

3 Measuring the sphaleron rate 

The evolution of the Chern-Simons number A'^cs over a time interval {to,t) can be defined 
using lattice electric and magnetic fields: 

6Ncs{t) = Ncsit) - Ncsito) = £2/' f d^^E'^B'^. (3.1) 

Unfortunately the topology on the lattice is not well defined, and using naive lattice scale 
E and B fields the right-hand side of equation (3.1) contains ultraviolet noise. This gives 
unphysical diffusion not connected with the sphaleron rate. The method of calibrated 
cooling [23] (see also [11]), offers a way out of the problem. It is based on the fact that 
at small enough lattice spacing sphalerons are large in lattice units, with a dominant 
length scale of order l/{g'^T).^ By applying a pre-determined amount of cooling (Langevin 
evolution without the noise) to the lattice gauge fields, the ultraviolet noise is eliminated, 
without compromising the long-distance topology of the configuration. At this point it 
is possible to evaluate the integral in (3.1) with only small errors. Cumulative residual 
errors are eliminated by periodically cooling all the way down to a vacuum configuration 
and correcting for deviation from integer values of SNqs between two vacua. This is 
schematically described in Figure 2. By adjusting the cooling parameters so that these 
deviations from integers always remain much smaller than unity, we also ensure that the 
cooling is sufficient to keep the measurement topological. 

''in ref. [12] the size of the sphaleron was estimated to be of order 5/(.g^r) in the pure SU(2) gauge 
theory. 
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Figure 2. Measurement of the Chern-Simons number evolution [23]. The solid circles show the 
configurations generated by the real-time evolution using the Langevin/heat-bath method. At fixed 
intervals, the configurations are cooled by the same amount in order to construct a cooled trajectory, 
where the UV noise is almost completely eliminated, allowing to calculate SNqs from (3.1). The 
cooling from vacuum to vacuum works as a test for residual errors: SNqs must then be close to an 
integer, the deviations from which are subtracted, thus avoiding the accumulation of errors. 



Cooling the original gauge fields close to the vacuum is computationally very expensive. 
The procedure is dramaticahy accelerated by blocking the lattice gauge fields by a factor of 
two after the UV noise has been sufficiently eliminated; this is repeated a couple of times 
until a minimum lattice size has been reached. For details, we refer to [23]. 

3.1 Sphaleron rate in the symmetric phase 

We calculate the sphaleron rate using two different, and complementary, methods. We start 
at high temperatures, above the cross-over, and go through the whole cross-over range into 
the "broken phase" by uniformly decreasing the temperature. At high temperatures we use 
standard canonical Monte Carlo sampling. As the potential barrier between consecutive 
Chern-Simons numbers is low, the probability distribution over Chern-Simons number is 
approximately flat. As we decrease the temperature, the sphaleron rate becomes exponen- 
tially suppressed and the canonical real-time method is too inefficient to resolve the slow 
rate. At this point we begin using multicanonical simulations. 

As an example, the evolution of the Chern-Simons number at rriH = 115 GeV and 
with T = 152 GeV (symmetric phase), 145 GeV (cross-over region) and 140 GeV (broken 
phase) is shown in Figure 3. In the symmetric phase the transitions are unsuppressed and 
it is straightforward to measure the diffusion rate. Around the cross-over temperature the 
probability distribution of AA'^cs becomes peaked around integer values and the transitions 
between these values become rapidly more suppressed. Finally, deep in the broken phase 
the rate goes down until we are not able to measure it with the real-time evolution method. 

3.2 Sphaleron rate in the broken phase: multicanonical method 

At low temperatures where the sphaleron rate is strongly suppressed, this can be measured 
using a multicanonical method similar to the one used in [23]. The calculation consists of 
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Figure 3. A^cs from a heat-bath trajectory (left), and the resufting probabiUty distribution (right), 
folded into the interval [0,1], at tuh = 115 GeV and T = 152 (top), 145 (middle) and 140 GeV 
(bottom). At high temperature, in the symmetric phase, the sphaleron transitions are unsuppressed, 
whereas at low T the transitions are so strongly suppressed that they do not happen in canonical 
simulations. 
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Figure 4. A heat-bath trajectory for Ncs, stih for mn = 115 GeV and T — 140 but now with 
multicanonical simulations (left), and the corresponding multicanonical probability distribution 
-Pmuca (right). 




Chern-Simons number Chern-Simons number N 



Figure 5. The physical distribution Pcan (left), after reweighting the result in figure 3.2 with the 
multicanonical weight function (right). 

two stages: 

i) the measurement of the potential barrier (probabilistic suppression) between two integer 
vacua, and 

ii) the calculation of the rate of the tunneling through the top of the potential barrier. 

Let us first look at the measurement of the potential barrier. This is the multicanon- 
ical stage of the computation. As is usually done in a multicanonical context, instead of 
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sampling configurations with the canonical weight 

Pcan(f/,<^)«e-^(^'*\ (3.2) 

we compensate for the strong suppression by adding a carefully-chosen weight function 
W{Ncs), which is a function of the order parameter, in this case the Chern-Simons number. 
The configurations {[/,<!>} are now sampled with the probability density 

Pmuca(C/,<^)oce-^(^'*)+^(^cs[^/]). (3.3) 
Defining the physical (canonical) probability distribution of the Chern-Simons number 

Pcan(iVcs) = / df/d<I> Pean[^7, m^hs " NcSp]), (3.4) 

it is clear that the corresponding multicanonical distribution is 

Pmuca(iVcs) =Pcan(iVcs)e^^'^°^^. (3.5) 

Thus, the probability suppression in multicanonical simulations vanishes if we choose 
W{N) = — lnpcan(-^) + const. This is not a particularly useful result because we do 
not know the canonical distribution a priori; indeed, that is the quantity we set out to 
compute with the multicanonical method. 

However, it is possible to calculate a good enough approximation for W by using an 
automatic iterative "self-learning" procedure. Here we follow the method presented in 
ref. [45], with the obvious modification of using a different order parameter. Essentially, 
during the learning stage the weight function is continuously modified in order to maximize 
the flatness of the total distribution of the Chern-Simons number. When the iteration has 
sufficiently converged, the resulting weight function W is then used in a production run. 

The physical (canonical) Chern-Simons probability distribution is now obtained from 
the measured multicanonical distribution using equation (3.5). An example of the mul- 
ticanonical evolution and the resulting probability distribution is shown in Figure 4, for 
niH = 115 GeV and T = 140 GeV; the same parameters as the lowest temperature in Fig- 
ure 3. We recall that in this figure the "evolution" cannot be interpreted as a physical 
evolution in real time. As we can observe, the distribution of the Chern-Simons number 
is now almost flat. The resulting canonical (physical) probability distribution is shown in 
Figure 5, together with the weight function W{Ncs) used. 

We obtain an estimate for the physical expectation value for a general observable A 
from multicanonical simulation from 

where the sums go over the configurations {U,^}i obtained from the simulation and on 
which the measurements are performed, Wi = W{Ncs[{U}i\) and Ai = A{{U, 

A couple of comments on the practical implementation of the multicanonical sampling 
using the Chern-Simons number are in order. Because now the Chern-Simons number 
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enters in the sampling weight for the configurations, it has to be evaluated in an unbiased 
manner; that is, the Chern-Simons number for a given gauge field configuration has to be a 
unique value (modulo 1), independent of the "history" of the configuration. The calibrated 
cooling described above has a small residual component which depends on the trajectory 
of the configurations. To ensure unbiased sampling, we do not use the cooled trajectory 
method but cool down to the vacuum at every measurement of A'^cs- More precisely, we first 
cool the gauge field a pre-determined amount (in most cases the cooling time is ~ 8.4 a^) 
in order to reduce the UV noise, and then start integrating E • B while cooling to vacuum. 
The frequent cooling now completely dominates the CPU-time budget of the simulation. 

The multicanonical probabilistic weight is implemented as an accept /reject step as 
follows: i) start with configuration A, with weight function Wa', ii) perform one heat- 
bath sweep through the lattice, producing provisional new configuration B; iii) measure 
Ncs{B) as described above, obtaining Wb- iv) The new configuration is accepted with the 
probability 



If the update is rejected, we start again at point i) with configuration A. The acceptance 
rate was around 50% at the lowest temperatures used, and increasing at higher tempera- 
tures. 

Obviously, the measured value of A'^cs depends on the amount of cooling applied before 
the measurement. Thus, the obtained probability distribution p(A'^cs) is also cooling depen- 
dent. However, this is completely cancelled by the dynamical rate measurement described 
in Section 3.3, so that the final rate is independent of the amount of cooling. Nevertheless, 
the right amount of cooling must be judiciously chosen for efficiency: insufficient cooling 
gives too noisy observables, whereas too much cooling takes one too far "downhill" from 
the original configuration towards the vacuum. In both cases the measured A'^cs is too far 
from the "true" value, and the result is that we do not observe a random walk for A'^cs- 
The situation becomes worse at large volumes and coarse lattices. 

3.3 Sphaleron rate in the broken phase: dynamical prefactor 

The multicanonical procedure described above gave us the probability distribution of the 
Chern-Simons number in the broken phase. We can now measure the tunneling rate fol- 
lowing refs. [23, 45]: 

1. Let us assume that we have done the multicanonical simulations and obtained the 
canonical (physical) probability distribution of the Chern-Simons number Pphys. ( A'^cs ) • 

2. We choose a narrow interval 1/2 — e/2 < A'cs ^ 1/2 + e/2 around the point that 
separates vacuum A^cs = from the vacuum A'^cs = 1- The relative probability of 
finding a configuration here is 




(3.7) 




(3.8) 
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1 ' I ' I ' I r 
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Figure 6. Two real-time trajectories starting from the same configuration. The final and initial 
configuration can either be the same (left) or different (right). The trajectories cross the central 
value of our order parameter A^cs = 1/2 several times, a fact we compensate for through the 
dynamical prcfactor (3.9). 



This is where we need multicanonical methods, as the probability of being on top of 
the barrier is extremely small, and to get a reliable estimate would take an imprac- 
tically long time with canonical sampling. 

3. Let us now take a random configuration from the canonical distribution but with the 
constraint 1/2 — e/2 < A'^cs < 1/2 + e/2; i.e. near the top of the potential barrier. 
Starting from this configuration, we now generate two real-time trajectories using the 
heat-bath dynamics, as described in Section 3.1. The trajectories are evolved until the 
Chern-Simons number falls near a vacuum value. Interpreting one of the trajectories 
as evolving backwards in time, we can glue the trajectories together at the starting 
point and obtain a vacuum-to-vacuum trajectory. The trajectory can either return 
to the starting vacuum or be a genuine tunneling trajectory, see Figure 6. Only the 
latter-type trajectories contribute to the sphaleron rate. 

4. We can obtain the tunneling rate by measuring lAA^^cs/^^l from the trajectories at 
the moment they cross the value A'^cs = 1/2. Here At is the time interval between 
successive measurements, and ANqs the change in Chern-Simons number. This 
characterizes the probability flux thorough the top of the barrier. We obtain the 
physical time difference from the relation between the heat-bath "time" and physical 
time, equation (2.12). 

5. If the tunneling trajectories would go straight across the top, the ingredients above 
would be sufficient to calculate the total rate. However, typically the trajectories 
"random walk" near the top of the barrier and can cross the value A'^cs = 1/2 several 
times. Because the trajectories were chosen starting from a set of configurations near 
the top of the barrier, this leads to overcounting: the evolution could be started at 
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any point the A^cs = 1/2 limit is crossed. This can be compensated by calculating a 
dynamical prefactor 



d = y , , (3.9) 

JVtraj ^ # crossmgs 

where the sum goes over the ensemble of trajectories, A'^traj is the number of trajec- 
tories, (5tunnci is if the trajectory does not lead to a change of the vacuum and 1 if 
it does, and crossings) is the number of times the trajectory crosses A'^cs = 1/2. 



With these ingredients, the sphaleron rate now becomes 



PdA^cs - 1/2| < e/2) 



AiVcs 



At 



d. (3.10) 



We note that the result is independent of e as long as e ^ 1. It is also independent 
of the frequency At with which the Chern-Simons number is measured: if we decrease 
the measurement interval, the trajectories become more jagged due to the random-walk 
nature of the heat-bath updates. This will increase the number of the crossings of the 
value A'^cs = 1/2 and hence decrease d. However, the latter is completely compensated 
by a corresponding increase in {\/S.Nq^/ lS.t\) . If the measurement interval At is small 
enough, random walk arguments imply d oc (At)^/^ and {\^Nq^/ /S.t\) oc (At)~^/^. This is 
corroborated by the numerical data. Thus, equation (3.10) has a well-defined continuum 
limit. 



4 Results 

We concentrate on two physically significant observables, the sphaleron rate and the Higgs 
field expectation value as functions of temperature at mn = 115 GeV and niH = 160 GeV. 
For both quantities we check for the finite volume and finite lattice spacing effects. 

4.1 Higgs field v{T) 

The gauge invariant Higgs condensate ((/>^(/)) is a direct probe of the phase transition or 
cross-over. At high temperatures it is close to zero, and at low temperatures it acquires 
an expectation value which grows as the temperature decreases. Because of the additive 
renormalisation, Eq. (2.8), the symmetric-phase value can become negative. The results 
obtained at /3g = 9, lattice size 32^, are shown in Fig. 7. Note that we display {(fP') in units 
of g^T"^, and goes to infinity as temperature goes to zero. 

At high temperature in the symmetric phase, (<I>^<I>) is close, but not quite identical, 
to zero. As the temperature is lowered, we enter the cross-over region where ($^$) grows 
rapidly. At smaller Higgs masses {rriH < 72 GeV), this rapid growth becomes a discontinu- 
ous jump, indicating a first-order transition [4]. Below the cross-over region, ($^<I>)/((7^T^) 
settles to an almost linear increase. 

For mu = 115 GeV we measured at selected temperature values while varying 
the lattice spacing by more than a factor of two (/3(5 = 6 . . . 16). The results are shown 
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Temperature (GeV) Temperature (GeV) 



Figure 7. The Higgs field expectation value for mn = 115 GeV (left) and 160 GeV (right). Squares 
are for /3a = 9, volume {L/aY = 32^. For niH = 115 GeV we have performed the continuum limit 
extrapolation at selected temperatures using a range of lattice spacings /3g = 4/(330) = 6 ... 12 
and extrapolating linearly to continuum. We observe that the /3g = 9 result deviates less than 8% 
from the continuum limit in the range of temperatures studied. The lines arc to guide the eye. 

in Figure 8. In this case we can reliably take the continuum limit by linear extrapolation. 
We observe that when (i;^^(/>) is small, the cut-off effects are very small, and at lowest 
temperatures T ~ 130 GeV the /3g = 9 result deviates from the continuum limit by less 
than 8%. The physical volume was kept fixed, at Lg^ ~ 14. We have checked that this is 
a large-enough volume so that the residual finite-volume effect is unobservable within our 
statistical accuracy. 

4.2 The sphaleron rate rdifr(r) 

The sphaleron rate F/T^ for Higgs masses 115 GeV and 160 GeV is shown in Fig. 9, using 
^(3 = 9 data. As expected, at high temperature in the symmetric phase, the rate becomes 
insensitive to the temperature apart from the trivial scaling. In this region the standard 
"canonical" real-time evolution is sufficient. As we proceed into the cross-over region there 
is a rather sharp turnover, with a drop of 10~^ over 5 GeV. In this region, both the canonical 
and multicanonical methods were used, and they agree very convincingly. Deeper in the 
broken phase, the decrease in the rate flattens out somewhat to a clean exponential drop- 
off, and, using the multicanonical approach, we were able to follow the rate over 10 orders of 
magnitude. For comparison, we have included the extrapolation from [26], expected to be 
valid deep in the broken phase. We see that the slope is correct, but that the central value 
of the rate is off by about an order of magnitude; or equivalently that the temperature axis 
is shifted by about 2 GeV for a Higgs mass of 115 GeV and about 5 GeV in the 160 GeV 
case. The data is shown for Pa = 9, where the majority of our simulations were done. 
Sphalerons are extended objects, and thus it is necessary to check the finite volume 
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Figure 8. Left: the Higgs field calculated for niH = 115 GeV and several temperatures, with 
decreasing lattice spacing a, but keeping the volume constant, according to Table 2. The black line 
is the continuum extrapolation. Right: an example of the continuum extrapolation at a single value 
of the temperature T = 143 GeV. 




Temperature (GeV) Temperature (GeV) 

Figure 9. The sphaleron rate for mn = 115 GeV (left) and 160 GeV (right) at f3a = 9. The shaded 
band is the theoretical estimate (plus error ranges extrapolated from lattice results in [23]) for the 
broken phase and the horizontal lines for the symmetric phase, as calculated in [26] . The canonical 
and multicanonical results agree within errors at high temperatures. 

effects. Using constant lattice spacing /3g = 9 and lattice sizes L/a = 16-54 we observed 
no systematic finite-size dependence within our statistical accuracy. Thus, we can be 
confident that L = 32 a ~ 14/^3 is safely large enough at all temperatures. This result 
is in agreement with ref. [12], where the volume dependence of the sphaleron rate became 
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negligible at L>5/g^ in pure SU(2) gauge theory. 

As we did with the Higgs field expectation value, we investigated the dependence of 
the sphaleron rate on the lattice spacing. We chose a set of six temperatures in the interval 
130 - 155 GeV and measured the rate at f^c = 4/ (fffa) = 6-16, while keeping the physical 
volume approximately constant: L ~ 3.5/30 = l^/^f. The lattice spacings and volumes are 
shown in Table 1, and the resulting sphaleron rates are shown in Figure 10. 

In the symmetric phase the lattice-spacing dependence is very mild. Deep in the broken 
phase the rate appears to decrease as a is decreased. This can be understood in the light 
of the increasing Higgs field expectation value at smaller a, see Figure 8. 

However, deep in the broken phase (T < 145 GeV) and for the smallest lattice spacings 
(/5g ^ 14) our multicanonical order parameter, cooled A'^cs (Figure 2), becomes ineffective 
and we are not able to obtain a sufficiently accurate measurement of the rate for the 
proper continuum limit. This is due to the increased noise in the measurement at smaller 
lattice spacings: when the amplitude of the noise is of order unity, a large fraction of the 
configurations with (measured) Nqs near half-integer value are actually some distance from 
the top of the tunneling barrier. Thus, only a small fraction of these configurations will 
lead to tunneling trajectories. 

The amount of noise can be reduced by applying more cooling before the measurement 
of Acs- However, cooling evolves the configuration towards one of the vacua (A'cs integer), 
and with too much cooling the measured order parameter does not track the true Chern- 
Simons number well enough for effective update. We emphasize that despite these issues 
the multicanonical method remains exact in the limit of infinite statistics; it is only the 
efficiency of the method which suffers. 

Because of this issue our statistics at small lattice spacing is severely restricted and 
we cannot obtain a reliable continuum limit. Thus, our final answer remains the /3g = 9 
result, where we have most of the data. However, what the data indicates is that the true 
continuum limit is probably a factor of 2-3 below the /3g = 9 result deep in the broken 
phase, which very likely makes the agreement with ref. [26] in Figure 9 better. 



a (i/gi) 


Pg 


L 


L/(3g 


0.67 


6 


20 


3.3 


0.5 


8 


28 


3.5 


0.44 


9 


32 


3.56 


0.4 


10 


36 


3.6 


0.29 


14 


48 


3.4 


0.25 


16 


56 


3.5 



Table 1. Lattice values for the continuum limit, a is the lattice spacing, /3g is defined in (2.4), 
L is the size of our volume. From the ratio L/ Pg we notice that we keep the physical proportions 
constant while we diminish the size of the lattice spacing a. 
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Figure 10. Left: the sphaleron rate calculated for mn = 115 GeV and with different lattice spacings 
a while keeping the physical volume constant. Right: an example of the rate at T = 143 GcV as a 
function of the lattice spacing. Here the continuum limit has been extrapolated by assuming that 
the dominant error of logF is linear in a. However, if we extrapolate T itself linearly in a the result 
is compatible with vanishing rate in the continuum limit. 




Figure 11. Left: the evolution of baryon number B (red) and lepton number L (black) in the pres- 
ence of a lepton number source, turned on at T = 200 GeV. The source is turned off at T = 170 GeV 
and 130 GeV respectively (blue dashed lines). Right: The evolution of the ratio (5.8), with the source 
turned off at T = 170 GeV (green) and 130 GeV (blue). Inserted: the evolution of the quantity 

vit). 

5 A sample Leptogenesis calculation 

To gauge the impact of using the correct sphaleron rate and Higgs expectation value 
through the transition, we solved the Leptogenesis equations (l.lO)-(l.ll), under some 
simplifying assumptions. 

As the sphaleron rate is well-known in the symmetric phase, we focus our interest on 
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the cross-over region, to investigate the efficiency of lepton-to-baryon number conversion 
through the newly-calculated sphaleron rate. In order to enhance the effect of the sphaleron 
rate suppression at cross-over temperatures, we study two limiting cases: one where the 
lepton- number source fi{t) was shut off well before the start of the cross-over, and one 
where we let the source active throughout. 

We took nc = 3 and assumed that all lepton species are equivalent Lj = L, i = 1, 2, 3, 
Li = 3L. The initial baryon and lepton numbers vanish L{t = 0) = B{t = 0) = 0. The 
source for the lepton number is therefore taken to be operational between T = 200 GeV 
and T = Tcut-oflf, 



Mt) = fit) 



fo 
2 



1 — tanh 



t cut-off 



T 



(5.1) 



2 GeV 

Since the equations are linear, the normalization of /o is arbitrary. This leaves only the 

source shut-off temperature Tcut off ^ free parameter, which we varied from 170 to 

130 GeV^, before and after the cross-over, respectively. 

The full expressions for 7(t) and r]{t) in (1.10-1.11) read [26] 



7(t)=n|p(^ 



1-X 



v{T) 



rdiff(r) 

2^3 ■ 



vit) 



(5.2) 



where T = T{t), and the functions 



pix) 



Xix) 



3 [65 + 136nG + 44n^ + (117 + 72nG)i 
2nG [30 + 62nG + 20n2, + (54 + 33nG)3 

4 [5 + 12nG + 4n| + (9 + 6nG)x'^] 
65 + 136nG + 44^2, + (117 + 72nG)x'^ ' 



(5.3) 



(5.4) 



We calculated the evolution of lepton and baryon number from temperature 200 GeV 
down to 130 GeV. Using that to a good approximation in the early Universe T cx 1/a, with 
a{t) the scale factor, we have that 



ZJ'T' 

dt ~ dT' 



(5.5) 



where H is the Hubble rate, given by the (radiation-dominated) Friedman equation^ 

= To . Q* = 106.75, Mpi = 2.43 x 10^^ GeV. (5.6) 



,2 _ ^^9*T^ 
90M2 ' 



Over the range of temperatures used here, H = (8.8 — 3.7) x lO"^'^. 

Once the source is turned off, and in the limit that L and B evolve much faster than 
V and ry, we can write the equations in terms of 1" = i? + •q'iL 



dlny _ 7(r) 
d\nT ~ H{T) 



(5.7) 



*We chose the 2 GeV width to mimic a fast shut-off of the source. 

^We ignored the effect of g* changing shghtly as the top quark begins acquiring its mass. 
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so that y = is enforced unless "y/H is too small. And so if the sphaleron rate is fast 
enough, we expect 

B(T) , , 

^ ' -1. (5.8) 



3L{T)rj{T) 



We say that the system is in "equilibrium" when this relation is obeyed. We note that 
r]{v{T)/T = 0) = 0.549... and r]{v{T)/T = oo) = 0.48. 

The evolution of B and L in time is shown in Fig. 11 (left). Starting from zero at 
T = 200 GeV, the introduction of the source leads to a growing L and, through sphaleron 
processes, growing B. This continues until the source is switched off; in the examples 
shown here Tcut-off = 170 GeV and Tcut-off = 130 GeV. For the early cut-off, both B and 
L level off to some asymptotic value. But even without switching off the lepton source, 
at a temperature around 143 GeV the sphaleron rate becomes inefficient, and the baryon 
number levels off. Lepton number is still sourced, but having no longer S as a sink, the 
growth of L becomes steeper. Tfreeze-out = 143 GeV corresponds to ^(t)/H ~ 10, and 
v{T)/T ~ 0.5. 

Fig. 11 (right) shows the "equilibrium" condition (5.8) in time. We see an initial 
transient, but from then on the system is nicely in equilibrium until the freeze-out at 
T = 143 GeV. In this case, where the source is still on (blue line), the equilibrium condition 
is obviously broken by the continued sourcing of L. Since L increases and B is constant, 
the ratio decreases. 

What is perhaps more surprising is the ratio (5.8) when there is no source. Because the 
v{t) increases through the transition, the equilibrium value r](t) also changes (see inset). 
A large enough sphaleron rate relative to the Hubble rate would adjust B relative to L 
to accommodate this evolving "equilibrium", but as is clear from Fig. 11, this does not 
happen. r]{t) just decreases and so the ratio becomes larger than 1. As a consequence, the 
asymptotic B and L obey 

B{x = oo) ~ r/(Tfrcczc-out)3L(rfrcczc-out) = 1-06 X 3L(rfrcczc-out)??(a; = oo). (5.9) 

Fig. 11 is based on the itih = 115 GeV data. We did a similar calculation for the 
niH = 160 rate, giving the same picture but with Tfreezo-out = 175 GeV. 



6 Conclusion 

In this paper we have presented the quantitative, non-perturbatively calculated sphaleron 
rate in the minimal Standard Model, using values of the Higgs mass still (marginally) 
allowed by experiment, 115 GeV and 160 GeV. At these Higgs masses, the electroweak 
transition is known to be an equilibrium cross-over. 

We first probed the temperatures in order to find the range at which the cross-over 
takes place. We performed lattice simulations in the symmetric phase with canonical Monte 
Carlo methods using a straightforward heat-bath update algorithm of the Higgs and gauge 
fields. Here we found that the rate is unsuppressed and follows a random walk in time. 
When lowering the temperature to reach the broken phase, the rate becomes suppressed 
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and in order to observe transitions between vacua we have to switch to multicanonical 
Monte Carlo and real-time simulations. 

We see a perfect match between the data points obtained with canonical and multi- 
canonical Monte Carlo, with a smooth overlap in the cross-over region. This gives us a 
cross-check on the validity of the results. And by varying the physical volume and the 
lattice spacing, we found that both the sphaleron rate and Higgs field are close to their 
infinite- volume and continuum limits. 

We obtained the sphaleron rate as a function of the temperature, showing an obvious 
similarity in the rate at the two Higgs masses mn = 115 GeV and 160 GeV, Figure 9. In 
both cases, the rates become exponentially suppressed in the broken phase with similar 
slopes, while having the same (~ 5 x 10^'' T^) asymptotic value in the symmetric phase. 
We point out that the cross-over range is clearly related to the Higgs mass: 130 — 147 GeV 
for rriH = 115 GeV, and 160 — 182 GeV for uih = 160 GeV. The cross-over range is clearly 
noticeable also in the plots for the Higgs field. 

Moreover, right at the beginning of the cross-over, the sphaleron rate drops at its 
fastest. Both in the plots for the sphaleron rate and for the Higgs field, the curves become 
steepest promptly after the cross-over kicks in, around 145 — 148 GeV for the smaller Higgs 
mass and 177 — 185 GeV for the bigger one. 

The behaviour of the sphaleron rate we found is in agreement with results quoted in 
the literature, in the range where they exist. Our asymptotic value is of the same order of 
magnitude as in [26]. Direct comparison in the broken phase shows that while the slope is 
the same, the value is off by an order of magnitude. 

Finally, we input the obtained sphaleron rate into a simple-minded model of Leptogen- 
esis. We found that although the rate cuts off exponentially at the transition, the freeze-out 
of the baryon and lepton number happens slightly later, when the rate is about 10 in units 
of the Hubble rate. For precision calculations of the generated baryon asymmetry in such 
models, it is therefore important to take the gradual shut-off of sphaleron processes into 
account. 
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